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^ ■ Abstract. We show that the KdV6 equation recently studied in [1,2] is equivalent to the 

Q ■ Rosochatius deformation of KdV equation with self-consistent sources (RD-KdVESCS) 
O . recently presented in [9]. The t-type bi-Hamiltonian formalism of KdV6 equation (RD- 
KdVESCS) is constructed by taking x as evolution parameter. Some new solutions of 
KdV6 equation, such as soliton, positon and negaton solution, are presented. 
Mathematics Subject Classifications (2000): 35Q53, 37K35. 
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1 Introduction 

I ■ 

^ ■ Recently, the 5 authors of [1] applied the Painleve analysis to the class of sixth-order 
OO . nonlinear wave equations, and found 4 cases that pass the Painleve test. Three of those 
^ ' cases correspond to previously known integrable equations, whereas the fourth one turns 
out to be new: 

{dl + 8u^d^, + AU^^) {Ut + + ^ul) = 0. (1) 

OO : 

O ■ This equation, as it stands, does not belong to any recognizable theory. In the variables 
V = Ux, w = Ut + Uxxx + 6ii^, ([I]) is converted to 

^ ■ vt + Vxxx + I2vvx - = 0, (2a) 

Wxxx + 8vwx + Awvx = 0, (2b) 

which is referred as KdV6 equation in [1] and regarded as a nonholonomic deformation of 
the KdV equation. The authors of [1] found Lax pair and an auto-Backlund transforma- 
tion for ([2]). They claimed that ([2]) is different from the KdV equation with self-consistent 
sources (KdVESCS) and reported that they were unable to find higher symmetries and 
asked if higher conserved densities and a Hamiltonian formalism exist for ([2]). 

In [2], Kupershmidt described ([2]) as a nonholonomic perturbations of bi-Hamiltonian 
systems. By rescaling v and t in ([2]), one gets 

Ut = Quux + Uxxx - Wx, (3a) 
Wxxx + 4:UWx + 2wux = 0, (3b) 
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which can be converted into a nonholonomic perturbations of bi-Hamiltonian systems [2] 



ut = B\^-^) - B\w) = B\^) - B\w\ (4a) 
B^{w) = 0, (4b) 

where 

B^ = d = a^, B^ = d^ + 2{ud + du) (5) 
are the two standard Hamiltonian operators of the KdV hierarchy, n = 2, and 

H^ = U. Ho = , H'i = — , • • • 

'23 2 

Then the author in [2] beheved that he could prove the integrabihty of KdV6 equation 
by constructing the infinite commuting hierarchy KdV„6 (HI) with a common infinite set 
of conserved densities. Some solutions for (j2j) were obtained in [1,3]. 

The soliton equations with self-consistent sources (SESCS) have attracted much at- 
tention (see [4]- [7]) and have important physical applications, for example, the KdV 
equation with self-consistent sources (KdVESCS) describes the interaction of long and 
short capillary-gravity waves [4] . The Rosochatius deformation of finite-dimensional inte- 
grable Hamiltonian system (FDIHS) also has important physical application, for example, 
the Garner-Rosochatius system can be used to solve the multicomponent coupled nonlin- 
ear Schrodinger equation [8]. We generalized the Rosochatius deformation from FDIHS 
to SESCS and presented many Rosochatius deformations of SESCS (RD-SESCS) in [9], 
such as RD-KdVESCS which stationary reduction gives rise to the well-known generalized 
Henon-Heiles system [10]. 

In this paper, we would like to answer the questions mentioned in [1]. We will first 
show that ([3]) is equivalent to the Rosochatius deformation of KdVESCS (RD-KdVESCS) 
presented in [9]. It is known [11, 12] that some soliton equations have both x— and t— 
type Hamiltonian formulation. However the Hamiltonian formulation for KdV6 equation 
(RD-KdVESCS) can not be written in usual way. We will formulate it as an infinite- 
dimensional integrable bi-Hamiltonian system with a t— type Hamiltonian operator by 
taking t as the 'spatial' variable and x as the evolution parameter as in the case of 
KdVESCS [13,14]. Since the KdV6 equation can be regarded as the KdV equation with 
no-homogeneous term and w is related to the square of eigenfunction, we may apply the 
method of variant of constant to find some new solutions of KdV6 equation starting from 
the known solutions of KdV equation. 

The present paper is organized as follows. We will first convert the KdV6 equation into 
RD-KdVESCS, and present extension of KdV6 equation in section 2. In section 3, we will 
describe RD-KdVESCS (KdV6 equation) and RD-mKdVESCS as a t— type Hamiltonian 
system by taking x as the evolution parameter, respectively. Then following the procedure 
given in [11]- [14] by means of the t— type Miura transformation relating these two 
Hamiltonian systems, we will construct the second t— type Hamiltonian structure for 
KdV6 equation (RD-KdVESCS) from the first Hamiltonian structure of RD-mKdVESCS, 



and present infinite chain of local commuting vector fields for KdV6 equation. Finally 
in section 4, starting from the solutions of KdV, we obtain many new solutions of KdV6 
equation, such as soliton, positon and negaton solution. 



2 KdV6 equation is equivalent to RD-KdVESCS 

By rescaling u and t and using the Galilean invariance of KdV equation, KdV6 equation 
([3]) can be rewritten as 

1 

ut = -{u^xx + 6uu^) - w^, (6a) 
Wxxx + 4(n - Ai)^:^. + 2wux = (6b) 

where Ai is a parameter. 

Set 

w = ^\ (7) 

then (6b) yields 

w^^^ + 4:{u - \i)w^ + 2wu^ = 2(p[ip.^^ + {u- Xi)(p]x + 6(Px[^xx + {u - Xi)(p] = 0, 
which immediately gives rise to 

(fixx + [u- Ai)v? = — , 

where /u is an integrable constant. So KdV6 equation is equivalent to 

Ut = ^{uxxx + 6uux) - {'^^)x, (8a) 
(Pxx + {u- Xi)lp = (8b) 

which is just the RD-KdVESCS presented in [9]. The Lax pair for (IHD reads [9] 

*^ -f(!M f/^f M M 



N-( ^+fW W 1 r9bl 



y_iLX-^-\ + ^^' if J 2 A - Ai V + # 
More generally, the multi-component extension of KdV6 equation is given by 

1 ^ 
Ut = -^{u^xx + Guu^) - ^ Wj^, (10a) 

i=i 

Wjxx + 4(m - Xj)wjr, + 2uxWj = 0, j = 1, 2, • • • , • • • A^. (10b) 



Under the transformation Wj = (fp fllUI) can be converted into the following RD-KdVESCS 



The Lax pair for ffTTj) is given by (9a) with 



1 ^ 
ut = -^{u^xx + Guu^) - (11a) 

i=i 

^jxx + {u- Xj)Lpj = j = 1, 2, ■ ■ ■ , A^. (lib) 



N 



9 2^ X - 



A' - f A - ^ - f + 1 E '/'I t I 2 ^ A - A, I V',^. + -V^.y^. 



(12) 



3 Bi-Hamiltonian structure of KdV6 equation 



In this section, we will follow the method in [11]- [14] to construct the bi-Hamiltonian 
formalism for KdV6 equation (RD-KdVESCS). First we will present the t— type Hamil- 
tonian formahsm for RD-KdVESCS and RD-mKdVESCS. For the RD-KdVESCS (jH]), 

set 



^u^x + - ip^ = c, qt = c^;, 
q = u, p= -]-u^, Q = <p, P = <Px, R= {Q,q,P,p,c) 



(13) 



then (IH]) becomes x— evolution equations and can be written as a t— type Hamiltonian 
system 



/ 



P 



ip 



(Ai-g)g + ^ 



\ 



Qt 



2^ 



= noV//i, 



(14a) 



,8H SH 5H SH SH, 



where V means variational derivative, VH — , , , , , 

oQ oq oP op be 

t — type Poisson operator Ho and conserved density Hi are given by 

/ 1 \ 

1 

Ho = -1 

-10 

\ 29t / 



, and the 



(14b) 



(14c) 



A 



The Rosochatius deformation of mKdV equation with self-consistent source (RD-mKdVESCS) 
is defined as [9] 



Gv'^Vx) + -(V9l(/?2)a 



(15a) 
(15b) 



Let 



q = v, p = ^Vx, Q = (fii, P = ^2, R= {Q,q,P,p,c)'^, 

then RD-mKdVESCS ( |T5|) can be written as a t— type Hamiltonian system 



(16) 



where VH 



( qQ + \iP \ 
2p 

f -QP- 2c 
V -Qt J 

5H 6H 6H 5H 6H 



(17a) 



^ - , ^_ , - , ^_ , -T^Y , and the t — type Poisson operator Hq 
oQ dq oP dp do 

and conserved density Hi are given by 
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The Miura map relating systems (IT^ to (|T71) . ie i? = M{R), is given by 

M: Q = Q, 

q = - f - 2p, 

P = XiP + qQ, 

p = \t-\c- \qp + \qp. 

c=Hi-qt = -ig' - 1^ + p2 + 2cq+ ^X^P' + qQP + 

which can be proved through direct calculations. 
Denote 

M' = 



(17b) 



(17c) 



(It 



where is the Jacobi matrix consisting of Frechet derivative of M, M'* denotes adjoint 
of M' . According to the standard procedure [11]- [14], applying the map M (fTSll to the first 
Hamiltonian structure of RD-mKdVESCS f|T7|) . we can generate the second Hamihonian 
structure of the RD-KdVESCS ([HD, 
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p 

-8p + 2dt 
(Ai-g)Q + ^ 

qdt + dtq J 



(19) 



and the bi-Hamiltonian structure for KdV6 equation or RD-KdVESCS (jHj) under the 
transformation ([7]) and (fT3|) is given by 



R-r = Hi 



8Hi 



ni- 



5H. 



c. 



(20) 



6R SR 

Since the t— type Possion operator Ho is invertible, we can immediately construct a 
recursion operator 

$ = ni(no)-^ 

lP9r' \ 

-Apd^^ + 1 
^[(Ai-Q)Q + ^]ar^ 

(21) 

which has the hereditary property [11]. Applying $ to the vector field Ki (IHl), we can 
generate the hierarchy of Hamiltonian commuting vector fields (symmetries) 
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and obtain a hierarchy of infinite-dimensional integrable bi-Hamiltonian systems 

R, = $"-1^1 = K^ = UoVHn = IliVHn-U 



(22) 



(23) 



for example 



( XiP + 2pQ + \qP 

2QP + qt 

-2pP + XfQ - UiQg + -Q^-cQ 



-\P'+Pt + \qQ'-\\iQ' 



\ 



4Q2 



2QQt + ct 

with the conserved functional densities H2 given by 



^7^-qpt+^cQ +-C . 



4 The new solutions of KdV6 

The KdV equation is 



ut = ^{uxxx + Quux). (24) 



The lax pair is 



ipxx + utp = Xijj, (25a) 

7/ II 

= -f ^ + (2 + ^)^- (25b) 

Let the functions (^1, 025''' ,0n be n different solutions of the system (|25!) with the 
corresponding A = Ai, A2, ■ ■ ■ , A„. We construct two Wronskian determinants from these 
functions: 

m = iy (01,- ■■ ,0!"^^\02,--- ,0^")), (26a) 

W2 = w{<Pu--- A''r\<p2,--- An,--- (26b) 

where rrii > are given numbers and 0^"'' := d'^(l)j{x, X)\x=Xy The generalized Darboux 
transformation of equation and system (1251) is given by [15] 

u = u + 2dllnWi, (27a) 



namely, system fl25|) is covariant with respect to the action of (1271) . For any initial solution 
of (IMl) . u and i/j are new solution of and fl25l) . Now we will take m = in what follows. 

4.1 soliton solution 

In f l26|) . let n = 1, mi = 0, A = Ai = and take 

01 (x, t, k) = coshQ, i/ji{x, t, k) = sinhQ, (28a) 
e = ^(x + ^kH) + a, Oi = y (x + ^klt) + a (28b) 

where a is an arbitrary constant. By using (1261) and (1271) . we obtain the single-soliton 
solution and the corresponding eigenfunction with k = ki for the KdV equation ([2] 



^2 

u = -^sech'^Gi, (29) 
f3ki 

tpi{x,t,ki) = —^sechQi, (30) 

where /S is an arbitrary constant as well. 

Since KdV6 equation can be considered to be KdV equation with non- 
homogeneous terms and w is related to the square of eigenfunction by ([7]), we may apply 



the method of variation of constant to find the solution of Eq.([n]) by using the solution u 
of Eg. (1211) and corresponding eigenfunction ipi. Taking a and /3 in (28b) and (150]) to be 
time-dependent functions a{t) and P{t) and using ([7]), and requiring that 



u 



-^sech^Qi, 



w = ipi{x, t, ki) = — - — seen Bi, 

©1 = y(x+^A;^t) + a(t), 

satisfy the Eq.Q. We find that a{t) can be an arbitrary function of t and 

4a'(t) 



Pity 

So the single-soliton solution of KdV6 equation 

kl 



ki 



(31a) 
(31b) 
(31c) 

(32) 



is given by 



u 



w 



-seek 9i, 



-a'{t)kiseeh'^Qi. 



(33a) 
(33b) 



Its shape is shown in figure 1. Notice that Gi contains an arbitrary t-function a{t). This 
implies that the insertion of sources into KdV equation may cause the variation of the 
speed of the soliton solution. So the dynamics of solution of KdV6 equation turns out to 
be much richer than that of solution of KdV equation. 
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(a) (b) 
Figure 1. The shape of single soliton solution for u and w when a{t) — 
~2t, fci = 1, t = 3. 

4.2 The first and second order of positon solution 



In 



set n = 1, mi = 1, A 



4 ■ 



Ai = — ^ and take 



01 (x, t, k) = sinQ, ipii^x, t, k) = eosQ, 
k 



e 



{x + Xi(k) 



—k'^t) ik — A;i)a, 



(34a) 
(34b) 



o 



where Xi{k) is a function that is analytic in the vicinity of the point k and has real Taylor 
expansion coefficients. By using (12^ , (ITTI) and we obtain first order of the one- 
positon solution and the corresponding eigenfunction with k = ki for the KdV equation 

m [15] 

— 16A;^sin0i(8sin9i + fci7cos6i) 
= (to2ei + A:i7)2 ' ^^^^^ 

^^("'^'^^^ = - 4.m29, + fe,7 ' ^''^^ 
ei = ^(x + xi(A;i)-^A;2t), (35c) 

7 = -8dkQ\k=ki = ^klt - 4(x + X2(fci)) + a, X2{ki) = [xi + 4:kdkXi{k)]k=ki (35d) 

where a, P are arbitrary constants. Similarly, by using ([7]) and the method of variation 
of constant we present first order of the one-positon solution for the KdV6 equation ([6]) 

— 16A;?sm0i(8sm6i + A;i7cos6i) , 
u = 7 ; — , 36a 

^ = 3klt-A{x + X2iki)) + a{t). (36c) 
(136|) implies that for fixed t and x ±oo, we have the asymptotic estimate 

u = -lsin2ei[l + 0{x-^)], (37a) 



w 



kia'it) 



sin^Qi[l + 0{x-% (37b) 



If X is fixed and t ±oo, the solution has the asymptotic behavior 

« = + )]' (38a) 

w = -^^stn'Q,[l + 0{t-% (38b) 

A positon solution as a function of x, u, and if have a second-order pole. This pole is 
situated at the point x = XQ{t) which oscillates around the point Xas{t) with the amplitude 
^, where Xasit) = —jklt — X2{ki) — The exact position of this pole can be determined 
by solving the following equation with S = kij, 

S = -4:sin^[6 - Akft + iki^Xi - X2) - kia{t)]. 
8 

So the positon solution of KdV6 equation ([6]) is long-range analogue of soliton and is 
slowly decreasing, oscillating solution. The shape and motion of the single positon is 
shown in figure 2. 
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Figure 2. The shape and motion of one-positon solution for u and w when a{t) = 
-2t, xi{ki) = 2ki, ki = 1. 

In order to find the second order of one-positon solution, we take © in (34a) to be 

e = ^(x + xi{k) - ^kH) -lik- kifa, 



(39) 



we have 



Wi^W{(t>i,dk(P,,df 



1 



128 



{-32sm^eicosei + kfYcosei + [Ukfu 



- 4^1 (4x + 4xi(A;i) + k^a - 2kix'l{ki))]sinQi}, (40a 
W2 = W{(l)i,dk(f>udlcf)i,jPi)\k=k, = -^kl{4sin2ei + ki^), 



hi 1 
©1 = —{x + xi{ki) - -kit), 7 



(40b) 

(40c) 
(40d) 



64 

ldkQ\k=k, = ^klt-A{x + X2{ki)), 

u = -4:dlQ\k=ki = 3/cit - 4:dkX2{k)\k=ki + «• 

By the similar process , by taking a = a{t) in (40d) we obtain the second-order of the 
one-positon solution from (7) and (27) 

u = 2dllnWu (41a) 

2a'{t) W2,2 



w — 



kf Vi ■ 



(41b) 



Its shape is shown in figure 3. 





rirv 


f\ 










-1 ) 










-2 










-31 










-4( 










-S( 










-6i 







-20 -10 



21. 



10 20 



(a) (b) 
Figure 3. The shape of second-order positon solution for u and w when 
a;i(fci) = 2kx, a{t) = -2t^, fci = 1, t = 2. 



1 n 



4.3 The first and second order of negaton solution 



In ( [26|) . set n = 1, mi = 1, A = ^, Ai = ^ and take 

01 (x, t, k) = sinhQ, ■0i(x, t, k) = coshQ, 

k 11 
e = -(x + xi{k) + -kh) + -(A; - ki)a. 
1 4 o 



(42a) 
(42b) 



We obtain the first order of one-negaton solution and tlie corresponding eigenfunction 
witli k = ki for the KdV equation ^ 



u 



— 16klsinhQi{8sinhQi — ki'ycoshQi) 



ipi{x,t, ki) 



(4sm/i26i - kijf 
AfiklsinhQi 
4sinh2Qi — ki'j ' 



h 1 
©1 = y(x + xi(A;i) + -k^t), 

7 = 84©|fc=fci = Skjt + 4(x + X2{ki)) + a, 
where a, (3 are arbitrary constants. 
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(43b) 

(43c) 
(43d) 
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Figure 4. The shape and motion of one-negaton solution for u and w when a{t) = 
-2i, xi{ki) = 2fci, ki = 1. 

Similarly, by using ([7]) and the method of variation of constant we present the 
first order of one-negaton solution for the KdV6 equation ([6]) 



1 1 



u 



w 



— 16klsinhQi{8sinhQi — kijcoshQi] 
(4sm/i26i - fci7)2 



(4sm/i2ei - k^y- 



7 = Sk^t + 4(x + X2iki)) + a{t). 

Similarly, negaton solution of ([6]) have second-order pole. The shape and motion 
negaton is shown in figure 4. 
Now we take G in (42a) to be 

k 11 

e = -(x + xi{k) + -kH) + -{k - kifa, 
z 4 8 



(44a) 

(44b) 
(44c) 
of the 



(45) 



we have 

W^ = W{(t)^,dk(t>i,dl4>i)\u=k, 



1 



128 



{'i2sinh'QiCoshQi - kfYcoshQi + [Ukfu 



- Aki{-4x - Axi{ki) + kia + 2kfx'l{ki))]sinhQi}, 

W2 = W^(0i, 5,01, 5^01, V^i)lfc=fci = ^kf{-Aszn2Q^ + k^^), 

o4 

ki 1 
Qi = —{x + xi{ki) + -kit), 7 = 8dkQ\k=ki 



3klt + A{x + X2iki)), 



V = 4dl<d\k=ki = 3kit + 4:dkX2{k)\k=ki + «. 
Similarly, we obtain the second-order of negaton solution from (7) and (27) 

u = 2dllnWi, 

2a'{t) W2.2 



w 



kf Wi' 



(46a) 
(46b) 

(46c) 
(46d) 

(47a) 
(47b) 



Its shape is shown in figure 5. 
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Figure 5. The shape of second-order negaton sohition for u and w when 
xi{ki) = 2fci, a{t) = 2t2, fci = 1, t = 5. 
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